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Abstract. 

Direct photon spectra are calculated within a transport+hydrodynamics hybrid approach, in 
which the high-density part of the transport evolution has been replaced by a 3+1-dimensional 
hydrodynamic calculation. We study the effects of changing the parameters of the two interfaces 
between the transport- and hydrodynamic descriptions on the resulting direct photon spectra. 



1. Introduction 

Heavy-ion physics is the principal tool to investigate the phase diagram of strongly interacting 
matter. Knowledge about the matter created in high-energy nuclear collisions is inferred by 
studying the remnants of the violent interaction. Since no more direct observation is possible, the 
experimental observations need to be complemented by theoretical calculations. Unfortunately, 
first-principle calculations of the whole collision are not feasible, so effective models need to be 
employed. 

Most remnants from heavy-ion collisions are themselves strongly interacting particles, which 
therefore are mostly emitted in the late stages of the fireball evolution. While they may carry 
information about the early stages collectively, such as flow patterns, the individual hadron 
carries only indirect information about the hot and dense phases. Direct photons are well-suited 
to probe the early stages of a heavy-ion reaction, since due to their small scattering cross- 
section they leave the reaction zone essentially unscathed. The experimental task of extracting 
direct photon spectra is challenging, because the decay of hadrons such as the 7r° create a 
huge background of decay photons. In this article, we investigate a model for direct photon 
emission [Ij with respect to the various parameters used therein. 

2. The Model 

We use the UrQMD v3.3 transport model [3J, Hj and its recent extension to substitute the 
intermediate high-density stage with an ideal 3+1-dimensional hydrodynamic calculation [5]. 
The early and late stages are still calculated in the transport model, which implements hadron 
and string degrees of freedom and uses PYTHIA [6] for scatterings at high momentum transfer. 
In the intermediate stage, hydrodynamic calculations are performed using three different 
Equations of State (EoS): A hadron gas EoS (HG-EoS [7]) is used, which has no phase transition, 
but the same the degrees of freedom as the transport description. A chiral EoS (x-EoS [8]) which 



has a cross-over phase transition to chirally restored and deconfined matter, and a MIT-bag 
model EoS (BM-EoS [9]) which has a strong first order phase transition with large latent heat 
to deconfined matter, are also used. 

The interfaces between transport and hydrodynamics are the main subject of the present 
investigations. Transport — > hydrodynamics: All particles and their momenta are transformed 
to baryon number-, energy- and momentum densities. In this process, the system is forced into 
local thermal equilibrium, because in ideal hydrodynamics only perfectly equilibrated matter 
can be described. The time at which this transition happens i s tart is, in standard setup, when 
the initial nuclei have passed through each other. This number depends on the (radii of the) 
incoming nuclei as well as the incident energy. We investigate the changing of this parameter 
from jio to 4to, with to being the standard value. At high energies, this time is too short to 
allow even for partial thermalization. The minimal time for switching to the hydrodynamic 
description has been set to t start = 0.6 fm, which is OAlto for Pb + Pfo-collisions at top SPS- 
energies (-E| a b = 158 AGeV) and more than 4to at ^/snn = 200 GeV. Therefore, the latter energy 
has been omitted from that investigation. 

Hydrodynamics —> transport: Besides the criterion when the transition from hydrodynamics 
back to transport happens, also the details of this transition can be changed. Usually, the 
transition happens when the system has diluted below e = 5eo, where Eq = 146 MeV/fm 3 
is the nuclear ground state energy density. This number is varied between £ cr j t = 2.5eo and 
£crit = 10. The scenario for the transition can either be that the criterion must be met by 
the part of the system that has the same z-coordinate (along the beam axis), whereafter this 
slice is transferred to the transport calculation. This is the default gradual scenario (for more 
information the reader is referred to jlOj). The other scenario, called isochronous scenario, 
requires the criterion to be met by the whole system, after which it is transferred to the transport 
calculation instantaneously. 

We will change the time for the first transition t sta n, the critical energy density for the 
second transition e cr jt and the scenario for the second transition (gradual vs. isochronous). 
During each of these changes, the other two parameters are kept at their default values. The 
calculations are done for central Pb + Pfr-collisions (6 < 5 fm) at E\ ab = 8, 45 and 158 AGeV 
and y/am = 200 AGeV. 

3. Results 

Figure Q] shows the invariant direct photon yield of calculations with varying start time of the 
hydrodynamic description. From left to right, the panels show the yield at 0.5 < p± < 1.5 GeV, 
1.5 < p± < 2.5 GeV and 2.5 <p± < 3.5 GeV. Each panel shows the yields for E\ ab = 158 AGeV, 
scaled by 10 (upper curves), E\ ab = 45 AGeV (middle curves) and E\ ab = 8 AGeV, scaled by 
1/10 (lower curves). At E\ ab = 8 AGeV and late switching times to > 2 fm, the whole system has 
an energy density below or close to the threshold for switching back to transport calculations 
e = 5eo, so that the hydrodynamic calculation only runs for a very short time. The yield of direct 
photons is maximal in calculations with t s tart ~ 0.75io for all systems. In calculations with large 
istart values, most of the high-density evolution of the system takes place in the transport phase. 
Here, chemical equilibration is reached only after some time, so that especially the irp — > jir- 
channel, which dominates the hadronic sources, is suppressed by the lack of /3-mesons in the 
unequilibrated cascade calculation with respect to the equilibrated hydrodynamic calculation. 

At the second transition from hydrodynamics to the transport calculation, the system remains 
in thermal equilibrium. The investigation shown in Figure [2] varies the critical energy density 
e cr it below which the system is described by transport calculations from e cr i t = 2eq, where the 
applicability of hydrodynamics seems very questionable to e cr j t = IOeq, where the use of a 
transport calculation might be questioned. We find a slight decrease of the photon spectra as 
the threshold is raised, consistent with the findings above. 
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Figure 1. (Color Online) Invariant direct photon yield of calculations with varying start time 
of the hydrodynamic description. The left panel shows the yield at 0.5 < p± < 1.5 GeV, the 
middel panel at 1.5 < p± < 2.5 GeV and the right panel shows the yield at 2.5 < p± < 3.5 GeV. 
Each panel shows the yields for E\ a b = 158 AGeV, scaled by 10 (upper curves), E\ a ^ = 45 AGeV 
(middle curves) and E\ 3 \> = 8 AGeV, scaled by 1/10 (lower curves) for HG-EoS calculations (blue 
solid lines), x~EoS calculations (orange dashed lines) and BM-EoS calculations (purple dotted 
lines). The scale is consistent within each panel. 
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Figure 2. (Color Online) Invariant direct photon yield of calculations with varying critical 
energy density for the transition from hydrodynamic description to transport description. For 
explanation of the panels and lines, see the description of Figure [TJ 



The influence of changing the transition scenario is investigated in Figure The yield in 
the gradual transition scenario is very similar to the yield in the isochronous transition scenario, 
albeit a small decrease is visible in most of the calculations. Since in the isochronous transition 
scenario the fraction of the system calculated with hydrodynamics is larger than in the gradual 
transition scenario, this is also consistent with the investigations on e cr i t and t s tart- 

4. Summary 

The parameter space of the UrQMD-hybrid model has been explored with respect to direct 
photon emission. The two major model parameters - the transition time from the initial 
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Figure 3. (Color Online) Invariant direct photon yield of calculations with varying the 
transitin scenario for the transition from hydrodynamic description to transport description. 
For explanation of the panels and lines, see the description of Figure [TJ 

transport calculation to the hydrodynamic calculation t s tart and the critical energy density for 
the transition time from hydrodynamics to the transport calculation e cr jt - have been varied 
over a range that includes all reasonable values. In all cases, we find that if a larger part of the 
evolution is calculated in hydrodynamics, the yield of direct photons increases. Only when the 
initial transition is done too early (small t s tart)j the direct photon yield is reduced. 

The increase of the photon spectra with longer hydrodynamic phase is expected in calculations 
with a Quark-Gluon-Plasma, which is only present in the hydrodynamic phase and emits more 
photons than the hadronic phase. In the hadronic phase, the increase is attributed to the higher 
number of p-mesons in the chemically equilibrated hydrodynamic calculation. 
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